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Abstract 

We consider a phenomenological extension of the minimal supersymmetric standard model 
which incorporates non-minimal chaotic inflation, driven by a quartic potential associated 
with the lightest right-handed sneutrino. Inflation is followed by a Peccei-Quinn phase tran- 
sition based on renormalizable supeipotential terms, which resolves the strong CP and fj. 
problems of the minimal supersymmetric standard model provided that one related parameter 
of the superpotential is somewhat small. Baryogenesis occurs via non-thermal leptogenesis, 
which is realized by the inflaton decay. Confronting our scenario with the current observa- 
tional data on the inflationary observables, the baryon assymetry of the universe, the gravitino 
limit on the reheating temperature and the upper bound on the light neutrino masses, we con- 
strain the effective Yukawa coupling involved in the decay of the inflaton to relatively small 
values and the inflaton mass to values lower than 10^^ GeV. 
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1 Introduction 

Recently non-minimal inflation (non-MI) [1], i.e. inflation arising in the presence of a non-minimal 
coupling between the inflaton field and the Ricci scalar curvature, TZ, has gained a fair amount of 
attention [2-8]. In particular, it is shown that non-minimal chaotic inflation based on a quartic potential 
[9] with a quadratic non-minimal coupling to gravity can be realized in both a non-supersymmetric 
[2, 3, 6] and a sypersymmetric (SUSY) framework [7, 8], provided that the inflaton couples strongly 
enough to IZ. In the latter case, the recently developed [8] superconformal approach to supergravity 
(SUGRA) greatly facilitates the relevant model building. In most of the models proposed, the inflaton 
is identified with the Higgs field(s) of the Standard Model (SM) or the next-to-MSSM {Minimal SUSY 
SM) [7, 8] - see also Ref. [10]. 

Motivated by the various attractive features of the MSSM [11] - such as the resolution of the 
hierarchy problem, the achievement of gauge coupling unification and the candidature of the lightest 
SUSY particle as cold dark matter - we consider it as the starting point of our investigation. Despite 
its successes, however, the MSSM fails to address a number of important issues. For instance, the 
strong CP and ji problems, the generation of the observed baryon asymmetry of the universe (BAU) 
and the existence of tiny but non-zero neutrino masses are some fundamental issues which remain 
open within the MSSM. For the resolution of these, it seems imperative to supplement the MSSM with 
additional superfields, which in the simplest cases are singlets under the SM gauge group, Gsm = 
SU{?,)c X S'C/(2)l X C/(l)y, so that gauge couphng unification is not disrupted. Consequently, new 
candidates (besides the Higgs boson) for driving non-MI arise. 

In Ref. [12-14] a resolution to the aforementioned problems of the MSSM was proposed within 
a framework that implements a Peccei-Quinn symmetry breaking phase transition (PQPT). In those 
models non-renormalizable superpotential terms are added, involving some singlets that develop vac- 
uum expectation values (VEVs) of the order of the PQ symmetry breaking scale. As a consequence, 
the // and the strong CP problems [15] of the MSSM can be simultaneously solved, and in addition 
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a new intermediate scale arises which generates Majorana masses for three right-handed (RH) neu- 
trinos, N?. The inclusion of N? is necessary so that the smallness of neutrino masses is explained 
through the well-known see-saw mechanism [16]. These same superfields can play an important role 
in the generation of the BAU via non-thermal leptogenesis [17-19]. This latter attractive possibiUty is 
invalidated, though, in the cases studied in Ref. [12, 13], where the PQPT follows a period of thermal 
inflation [20] that leads to a very low reheating temperature. An enormous entropy production occurs, 
diluting any preexisting, non-thermally created, lepton asymmetry. This dilution can be avoided if we 
adopt the scheme of Ref. [14] but then, the PQ field cannot be zero during inflation - see below. 

On the other hand, non-thermal leptogenesis can be enhanced if the scalar component, N, of the 
lightest is the inflaton itself as firstly proposed in Ref. [21]. In this case the branching ratio of the 
inflaton decay (which now triggers leptogenesis) into a lepton plus a Higgs boson is [19] maximized. 
However, iV-inflation, in its simplest realization [13,21,22], is of the chaotic type - for other scenarios 
see Ref. [23,24] - and therefore, trans-Planckian inflaton-field values are typically required to allow for 
a sufficiently long period of inflation. The implementation of inflation then necessitates the adoption 
of special types of Kahler potential, as in Ref. [18], so that SUGRA corrections are kept under control 

- for other proposals related to chaotic inflation with a quadratic potential, see Ref. [25]. Moreover, 
minimal chaotic inflation driven by a quartic potential seems [26] to be ruled out by the fitting to the 
seven-year data of the Wilkinson Microwave Anisotropy Probe Satellite (WMAP7), baryon-acoustic- 
oscillations (BAO) and Hubble constant (Hq) data [27]. 

In this paper we construct a model of non-minimal N inflation (non-MiVl ) retaining the successful 
ingredients of the picture above. To this aim, N (the lightest RH sneutrino) is coupled to one of the 
PQ fields, which can be confined to zero during inflation - see Ref. [8]. We then show that the model 
naturally leads to non-MI within SUGRA, provided that a particular parameter of the superpotential 
is sufficiently small. Sub-Planckian values of the inflaton field are allowed in a wide range of the 
parameter space, and the adopted type of Kahler potential is more or less well-motivated. Also the 
inflationary observables turn out to he within the current data. The non-MiVI is followed by a PQPT 
driven by renormaUzable superpotential terms as in Ref. [28, 29], whereas the p, parameter of the 
MSSM can be generated from the PQ scale via a non-renormalizable term as in Ref. [12-14]. The 
reheating temperature is determined exclusively by the decay of N and is high enough (> 100 GeV) 
so that non-perturbative electroweak sphalerons are operative and, consequently, non- thermal [17] 
leptogenesis and the subsequent generation of the BAU can be reaUzed. As usually in similar models 

- cf. Ref. [19, 22-25, 30] - consistency with the constraint on the gravitino (G) abundance [31-33] 
requires a relatively small effective Yukawa coupling constant (lO~® — 10~^). The smallness of this 
coupling though may be explained through a broken flavor symmetry [19, 30]. 

Below, we present the basic ingredients of our model (Sec. 2) and describe the inflationary (Sec. 3) 
and post-inflationary dynamics (Sec. 4). We then restrict the parameters of our model (Sec. 5) and 
summarize our conclusions (Sec. 6). Details concerning the formulation of non-minimally coupled 
scalar fields within SUGRA are presented in the Appendix. Throughout the text, we use natural units 
for Planck's constant, Boltzmann's constant and the speed of light {h = c = = 1); the subscript 
of type , X denotes derivation with respect to (w.r.t.) the field x (e-g-, ,xx — (^"^ /d'}^)', charge conju- 
gation is denoted by a star and log [In] stands for logarithm with basis 10 [e]. Finally, we follow the 
conventions of Ref. [34] for the quantities related to the gravitational sector of our model. 

2 IVIODEL Description 

We focus on a PQ invariant extension of the MSSM, inspired by Ref. [12, 13], which links the gener- 
ation of intermediate masses for RH neutrinos, Nf, with a PQPT. Besides the (color) anomalous PQ 
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Table 1: The representations under Gsm and the extra global charges of the superhelds of our model. 



symmetry U{1)pq, the model possesses also an anomalous R symmetry U{1)r, and the baryon num- 
ber symmetry U{1)b- The PQ symmetry U{1)pq is spontaneously broken at the PQ breaking scale 
fa ~ (10^° - 10^2) GeV (which coincides with the axion decay constant - for a review see Ref. [35]) 
via the VEVs acquired by two Gsm singlet left-handed superfields X and X. The representations 
under Gsm and the charges under the global symmetries of the various matter and Higgs superfields 
are listed in Table 1 . 

In particular, the superpotential, W, of our model naturally sphts into two parts: 

W = Wmssm + Wnpq, (2.1) 

where VFmssm is the part of W which contains the usual terms - except for the n term - of the MSSM, 
supplemented by Yukawa interactions among the left-handed leptons and N?: 

W^MSSM = hEijeiLjHd + hDijdlQjHd + huijulQjHu + hNijN^LjHu- (2.2) 

Here, the group indices have been suppressed; the i-th generation SU{2)i^ doublet left-handed quark 
and lepton superfields are denoted by Qi and Lj respectively, whereas the S'[/(2)l singlet antiquark 
[antilepton] superfields by ti? and di^ [e? and Nf] respectively. The electroweak Higgs superfields, 
which couple to the up [down] quark superfields, are denoted by Hu [H^]. 

On the other hand, VFnpq is the part of W which is relevant for non-MA^I, the generation of the 
Majorana masses for Nf, the spontaneous breaking of U(l)pQ and the generation of the /x term of the 
MSSM. It takes the form 

T^NPQ = XiXNfNt + XaP{XX - fl/4.) + A^^^^^, (2.3) 

where mp ~ 2.44 • 10^^ GeV is the reduced Planck scale and P is a Gsm singlet left-handed superfield 
involved in the breaking of C/(1)pq. The parameters \a and /„ are made positive by field redefinitions. 
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Moreover, we chose a basis in the Ni — Nj space where the coupling constant matrix A is real and diag- 
onal. In order to produce the CP-violation necessary for leptogenesis, we include three N?. Assuming 
that are strongly hierarchical, i.e. Ai = A ^ A2, A3, the scalar components of the two heavier 
roll to their minima fairly quickly, since their potential is steeper - especially if we further assume that 
these are minimally coupled to gravity in contrast to the lightest one. Thus the scalar component, N, 
of the lightest of the Nf's controls the relevant slow -roll dynamics and it can therefore be identified as 
the inflaton. On the other hand, the outcome of leptogenesis is governed by the one from N[ with the 
smallest decay rate. In the following, we concentrate on the case where N drives both non-MTVI and 
leptogenesis. Therefore, the three generation model can be simplified to an effective one-generation 
model in the Nf sector, with the only remnant of the other two generations being a non-vanishing 
CP-asymmetry for leptogenesis - see Ref. [21,23]. Henceforth we suppress family indices. 

According to our general discussion in the Appendix - see Eq. (A. 10) - the implementation of non- 
MNl within SUGRA requires the adoption of a frame function, $7 j^^, related to the Kahler potential, 
Kj^^, of the following form 

= -3e-^^./3-^ = _3 + ^ + ^ + ^ + ^-^^^-^(iV2 + N*A , (2.4) 

rup rup THp rrip rrip Amf, \ / 

where the complex scalar components of the superfields P, X and X are denoted by the same sym- 
bol and the coefficients kx and /c^ are taken, for simplicity, real. Comparing this expression with 
Eq. (A.9), we remark that we adopt the standard non-minimal coupling for the inflaton, N, i.e. F = 
kj^N'^/Anip, and we added the sixth term in the RH side (RHS) in order to cure the tachyonic mass 
problem encountered in similar models [7, 8] - see Sec. 3.1. Note that F breaks explicitly the imposed 
R and PQ symmetries and causes a dependence of O on the phase 6 of N which can be written as 
N = ae /V2. As we show in Sec. 3.1, the model admits stable inflationary trajectories along which 
6 is stabihzed at = 0. When ~ 0, the choice kj^ > 1 ensures the positivity of the scale function, 
— J7^j/3, even for relatively large values of a. 

In the limit where mp tends to infinity, the matter sector decouples from gravity, and we can obtain 
the SUSY Umit, FsuSY, of the SUGRA potential, V^j. This turns out to be 

^ \ ^ _ 2 

v2|A7'|2| \2|i:)|2A|-y"|2 " ~ 



V^susY = [4.X'\N\' + Xi \P\'j \X\' + XN' + A„PX + A^ \XX - f'jA\\ (2.5) 

From the potential in Eq. (2.5), we find that the SUSY vacuum lies at 

(TV) =0, (P) ~ and | | = 2| (X) | = 2| (X) | = /„, (2.6) 

where we have introduced the canonically normalized scalar field (j)x = 2X = 2X. Note that, 
since the sum of the arguments of {X), {X) must be 0, X and X can be brought to the real axis 
by an appropriate PQ transformation. Moreover, after including soft SUSY breaking terms, (P) can 
become [28] of order 1 TeV and the minimization of Vp in Eq. (2.5) requires that {X) = {X). Needless 
to say that for field values greater than /„, the TeV-scale soft SUSY breaking terms can be safely 
ignored during the cosmological evolution. After the spontaneous breaking of [7pq(1), the first term 
of the RH side (RHS) of Eq. (2.3) generates intermediate scale masses for the A^f s Mi ~ X-ifa and, 
thus, seesaw masses [16] for the light neutrinos, whereas the third of the RHS of Eq. (2.3) leads 
to the fi term of the MSSM, with ~ A^ 1(^)1^ /'mp. Both scales are of the right magnitude if 
\{X)\ = fa/2 ~ 5 • 10^1 GeV, Ai ~ 1 and A^ ~ (0.001 - 0.01). 

In conclusion, Wnpq leads to a spontaneous breaking of C/pq(1). The same superpotential Wnpq 
also gives rise to a stage of non-MTVI and a PQPT, as analyzed in Sec. 3. An indication for such a 
possibiUty can be seen by examining Vsusy in Eq. (2.5), which becomes 

VsvsY = X^N'^ + Xlf^/ A along the direction X = X = 0. (2.7) 
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Clearly, for N S> fa, Vsusy tends to a quartic potential. Therefore, Wnpq can be employed in 
conjunction with Kj^^ in Eq. (2.4) for the realization of non-MA^I along the lines of Ref. [8]. Moreover, 
for lower iV's, Vsusy takes an almost constant value, which can drive a PQPT. 

It should be mentioned that the non-minimal gravitational coupling, instanton and soft SUSY 
breaking effects explicitly break U{1)r x ?7(1)pq to a discrete subgroup. It is then important to 
ensure that this subgroup is not spontaneously broken by {X) and (X), since otherwise cosmologi- 
cally disastrous domain walls are produced [36] during the PQPT. Note that U {1)r x U{1)pq is also 
broken during non-MA^I due to the non-zero N, but it is restored in the SUSY vacuum. The explicitly 
unbroken subgroup of C/(1)r x C/(1)pq can be deduced from the solutions of the system 

4r + 2p = (mod 2tt) , 4r = (mod 2tt) and - 12(r + p) = (mod 27r) , (2.8) 

where r and p are the phases of a [/ (1)^ and U{1)pq rotation respectively. Here we took into account 
that: (a) the R [PQ] charge of the langrangian term caused by the non-minimal gravitational coupling 
is 4 [2]; (b) the R charge of W and, thus, of all the soft SUSY breaking terms, is 4 and (c) the sum 
of the R [PQ] charges of the SU{3)c triplets and antitriplets is —12 [—12]. We conclude, therefore, 
that the explicitly unbroken subgroup is Z4 x Z2. It is then easy to check that this subgroup is not 
spontaneously broken by {X) and {X), since the relevant condition 

2p = (mod 27r) , (2.9) 

is satisfied automatically as a result of Eq. (2.8). Consequently, cosmologically disastrous domain 
walls are not produced during the PQPT and so, there is no need to further extend [37] the particle 
content of the model - cf. Ref. [29]. 

3 The Inflationary Epoch 

3.1 Structure of the Inflationary Action 

Inserting Eq. (2.4) into Eq. (A. 8), we can write the action of our model in the JF (Jordan frame) as 
follows 

= J d^xV^ (^^p%i^ + S^-pd^rd'^ct^*'^ - Q.^^A^A'^/ml - V^^ , (3.1) 

with 0" = N,P,X and X. Also V^j = ^|^jV^i/9 where V^j is the EE (Einstein frame) E-term 
SUGRA scalar potential, which can be obtained from Wnpq in Eq. (2.3) - without the last term of the 
RHS - and in Eq. (2.4) by applying Eq. (A.3). Along the direction P = X = X = 0, V^j and 
= — 3/ take the forms 

^ATio = ™P with / = 1 - + ( - + C7^ j xl cos 29 and c^^ = -- + -f- (3.2) 

Here MpQ = /a/2mp and Xa = a/mp. Recall also that we set N = ae^^ /■\/2. From Eq. (3.2), we 
can easily verify - see also the small fluctuations analysis below - that for given a, 6 = (modulo 
vr) minimizes V^j. For 9 = and cjz ^ I, S in Eq. (3.1) takes a form suitable for the realization 
of non-MA^l. Then we can set = 0, and more importantly V^j develops a plateau - note that 
Mpq <C 1. The constant potential energy density V^jq and the corresponding Hubble parameter H^j^ 
along the trajectory for which non-MA^I can take place are given by 
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In order to check the stability of the direction P = X = X = 6 = w.r.t. the fluctuations of the 
fields 6,P,X and X, we expand the latter three in real and imaginary parts as follows 

„ Xl+iX2 ^ Xi+iX2 Pl + iP2 

Performing a Weyl transformation as described in Eq. (A.4), we obtain [38] 

+^9^''129f^xd.X-Vj^^ , (3.5) 

with x = xi,X2, xi,X2,pi,P2, and we also take into account that <^ f^^ and ffi <^ f^a for ^ ~ 0. 
Note that we keep only terms up to quadratic order in the fluctuations 0, x and their derivatives in 
Eq. (3.5). Along the trajectory P = X = X = 9 = 0, f becomes a function of a, and so we can 
introduce the EE canonically normalized fields, a, 9 and as follows [5,8] 

{§h''=rl"'KW' -oj-Tf- 0.6) 

Taking into account the approximate expressions for a - where the dot denotes derivation w.r.t. the 
cosmic time t - J and the slow-roll parameters rj, which are displayed in Sec. 3.2, we can verify 

that, during a stage of slow -roll non-MA^I, 6 ~ Ja9 since J a ~ Vdmp, and x — x/V7- For 
the latter, the quantity f/f^^^, involved in relating x to x, turns out to be negligibly small, since 
///3/2 = f^^a/f/'^ = -Xy^\mp/2V3cTz. Therefore the action in Eq. (3.5) takes the form 



J d^x^ i^\mln + ^^'^^ ^ d^^d,4> - % j , 



where (j) stands for cr, 9, xi,X2,xi,X2,pi and p2. 

Along the inflationary path, we can easily check that the first derivatives of V^j w.r.t. cf) are equal 

to zero. The curvature of V^j w.r.t. 9 can be studied separately since d'^V^^/d9dx = 0. Moreover 
the stability of the path P = X = X = 9 = {) w.r.t. the fluctuations of is automatic, since the mass 
squared of 0, m\, turns out to be positive. Indeed we find 

, \^ml{l + &cn)xi ^ X^ml _ 

Qj2p - 3^ -^^ATIO- ^^-^^ 

The two 3x3 mass squared matrices M\ = {ci^Vf^-^/ dxadXjS^ ^ with A = 1 and Xa = Pi,xi,xi or 
A = 2 and Xa = P2,X2,X2, have the following eigenvalues 

„| = (l^±4Zi|%Z^ ^ (3.9) 

corresponding to eigenstates xi (or X2) and yi± = (pi ± xi) / \/2 (or y2± = (p2 ± X2) / \/2) respec- 
tively. The considered inflationary trajectory, P = X = X = 9 = 0, isa stable valley of local minima, 
provided that m| > and ■mg_^ > 0, i.e. 



Tn / 3A 

(a) cr > (Tie = ^ and (b) a > a2c = mpJ - — with Xa < X/3cn- (3.10) 

Vokxcn V A - 6XaCTz 
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In practice the condition of Eq. (3.10b) is much more restrictive than Eq. (3.10a) for kx ~ 1- Indeed, 
from Eq. (3.9), it is evident that kx ^ I assists us to achieve m| > - in accordance with the results 
of Ref. [8]. On the other hand, given that for a < mp, we need cji S> 1, Eq. (3.10b) requires a clear 
hierarchy between A and Aa, e.g. for c-ji c±. 10^, we need Aa/A < 10~^. This ratio can be slightly 
increased (almost one order of magnitude) if we include additional terms such as /cji^^|iVp|Xp or 
kpff\N\'^\P\'^ in the Kahler potential in Eq. (2.4). Since the resulting increase of Aa has no significant 
impact on our results, we choose to stick with the most minimal possible Kahler potential needed 
for the viability of our model, avoiding more complications. We have also numerically verified that 

> ifjyp >H^ and mig^ >H^, during the last 50 — 60 e-foldings of non-MTVI, and so any 
inflationary perturbations of the fields 9, xi^2 and yi.2± are safely eliminated. 

The constant tree-level potential energy density in Eq. (3.3) causes SUSY breaking, leading to the 
generation of one-loop radiative corrections, which can be calculated by employing the well-known 
Coleman- Weinberg formula [39]. We find 

= ^ [rnpn ^ + 2m|ln ^ + 2mi^ In !^ + 2mi_ In !^ - Am' In ^ j (3.11) 

where A is a renormahzation mass scale and m = \/2\m,pXfj/ p^'^ is the eigenvalue of the fermion 
matrices. As we verified numerically, V^c has no significant effect on the inflationary dynamics. This 
is because the slope of the inflationary path is generated at the classical level - see the expressions for 
?and ?7 below - and so, the contribution of to V^j remains subdominant. 

Based on the action of Eq. (3.7) with V^j ~ ^mo + we can proceed to the analysis of non- 
MiVl in the EF, using the standard slow-roll approximation [40,41]. It can be shown [42] that the 
results calculated this way are the same as if we had calculated them using the non-minimally coupled 
scalar field in the JF. 

3.2 The Inflationary Observables 

According to our analysis above, when Eq. (3.10) is satisfied, the universe undergoes a period of slow- 
roll non-MATI , which is determined by the condition - see e.g. Ref. [40,41]: 

max{e'(cr), |r/((j)|} < 1, where 

and ?? = mf. ^' =^\^ ^ — l_ ~ _^ — £_. (3.12b) 

Here we employ Eq. (3.3) and the following approximate relations: 

J~V6mp/a, %^^~AV|/4a3 and -Z\^ml/c^a\ (3.13) 

The numerical computation reveals that non-MA^"! terminates due to the violation of the ? criterion at a 
value of a equal to erf, which is calculated to be 

e{ai) = l => ai = {A/3fmp/^- (3.14) 

Note, in passing, that for a > the evolution of a - or cr via Eq. (3.6) - is governed by the equation 
of motion 

3H^^& = -i/^j ^ ^H^^J^a = ^ & = -Xml/sVSc^a. (3.15) 
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Using Eqs. (3.12a), (3.12b) and (3.15), we can derive the expression for f / f^^^ given above Eq. (3.7). 

The number of e-foldings, A^*, that the scale fe* = 0.002/Mpc suffers during non-MA'^1 can be 
calculated through the relation 

N,= \ dd^ = \ da J2^, (3.16) 

where cr* [a*] is the value of a [a] when A;* crosses the inflationary horizon. Given that erf ^ a*, we 
can write cr* as a function of iV* as follows 

The power spectrum Pji of the curvature perturbations generated by a at the pivot scale is 
estimated as follows 

pi/2 ^ 1 %(^*)^/^ ^ l^(^.)l ^ Aa,^ ^ V2A^ ^ 

^ 2x/37rm3 2^/3 ^mf, 8^/27rm2 127rc7^ ' 

where Eq. (3.17) is employed to derive the last equality of the relation above. At the same pivot scale, 
we can also calulate the (scalar) spectral index, rig, its running, Og, and the scalar-to-tensor ratio, r, via 
the relations: 

ns = 1-6?* + 25?, ~ 1 - 2/iV*, (3.19a) 
«s = ^ (4??2 - (ns - 1)') - 2^* ^ -2/iV2 (3.19b) 
and r = 16e* ~ 12/N^, (3.19c) 



where ^ = ?np^i s^^ni sffj/^i^i ~ m^v'Zerj^fj/ J + 2^ and the variables with subscript * are 
evaluated at a = cr*. Comparing the results of this section with the observationally favored values, we 
constrain the parameters of our model in Sec. 5. 



4 The Post-Inflationary Evolution 

A complete SUSY inflationary scenario should specify the transition to the radiation dominated era 
and also explain the origin of the observed BAU consistently with the G constraint. These goals can 
be accomplished within our set-up, as we describe in this section. The basic features of the post- 
inflationary era of our model are exhibited in Sec. 4.1. A more precise analysis of the evolution during 
this era can be obtained by solving numerically the relevant Boltzmann equations, as in Sec. 4.2. 
Finally, useful analytical expressions reproducing accurately our results are presented in Sec. 4.3. 



4.1 The General Set-up 

When non-MiVl is over, N undergoes a very short period of fast-roll until it reaches its critical value 
C2c/V2- Afterwards, N and the PQ system (comprised by X, X and P) fall to their SUSY minimum 
values, acquiring masses and mpQ respectively, which can be computed from Vsusy in Eq. (2.5). 
These are given by 

(a) = A/a and (b) mpQ = \afa/V2- (4.1) 

Note that due to hierarchy between A and Aq established in Eq. (3.10b), m^^ > mpq. Consequently, the 
post-inflationary energy density of the universe is dominated by the N condensate which undergoes a 
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phase of damped oscillations about the SUSY vacuum, when i7 ~ m^, and decays [21] predominantly 
into Hu + Lor H* + L*, via the tree-level couplings derived from the last term in the RHS of Eq. (2.2). 
The initial energy density of this oscillatory system is estimated by pii ~ 3mpm?~, corresponding to 

Hii ^ rrij^. The decay temperature of N, T^, which coincides with the reheating temperature, Trh, in 
our model - see below - is [44] given by 



= CT Y^r^mp with CT = ( 



72 



(4.2) 



where counts the effective number of relativistic degrees of freedom at temperature T^. We find 
5* ~ 240 for the MSSM spectrum plus the particle content of the superfields P, X and X. Also is 
the decay width of N given by 

r^ = ^^eff"ijV "^^^^^ ^eff = \/^i\hNuf (4.3) 

is an effective Yukawa coupling, linked to the Ught neutrino masses, and can be considered as a free 
parameter. 

The aforementioned two channels for the N decay have different branching ratios when CP conser- 
vation is violated. Interference between tree-level and one-loop diagrams generates a lepton-number 
asymmetry [19, 34] which, for a normal hierarchical mass spectrum of light neutrinos, reads 

3 m,y3 

Here |5eflf| < 1, which is treated as a free parameter in our approach, represents the magnitude of CP 
violation; 1711,3 is the heaviest neutrino mass and we take (Hy) = 174 GeV (adopting the large tan^ 
regime). If Trh < rrij^, the out-of-equilibrium condition [19] for the implementation of leptogenesis 
is automatically satisfied. The resulting lepton-number asymmetry after reheating can be partially 
converted through sphaleron effects into baryon-number asymmetry. However, the required Trh must 
be compatible with constraints for the G abundance, Yq, at the onset of nucleosynthesis. 

On the other hand, the system consisting of the two complex scalar fields P and {SX + SX) / \/2 
(where SX = X — fa/2 and 5X = X — fa/^,) enters into an oscillatory phase about the PQ minimum 
and eventually decays, via the non-renormalizable coupling in the RHS of Eq. (2.3), to Higgses and 
Higgsinos with a common decay width Fpg [41] and a corresponding decay temperature Tpg given 

by 

TpQ = CT^/TpQmp where Tpq = (2^) ^^'^^ 

Note that due to the hierarchy between A and A^ in Eq. (3.10), the decay of X to ATP's is kinematically 
forbidden. Due to the same fact, Tpq turns out to be quite suppressed, and so a possible domination 
of the PQ oscillatory system could dilute any preexisting Yl and Fg. However, {X) <C mp - in 
contrast to the VEVs of moduli occurring in superstring theory [20] which are of the order of mp - 
and therefore, the initial energy density of the X and X's oscillations, p2i, is reduced w.rt. the energy 
density of the universe at the onset of these oscillations, pijpq - cf. Ref. [19, 30]. Namely 

P2i c± mpQ|(X)|^ < pHPQ = 3m|.m|.Q for H ~ mpq. (4.6) 

This is a crucial point since it assists us to avoid any dilution of the produced Yl (and Yq) at T = T^, 
as we show in the following. 



4 The Post-IaHationary Evolution 
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4.2 The Relevant Boltzmann Equations 

The energy density, pi [p2], of the oscillatory system with decay width [rpg], the energy density 
of produced radiation, pR, the number density of lepton asymmetry, ul, and the one of G, uq, satisfy 
the following Boltzmann equations -cf. Ref. [21,29,33,34,43]: 

/9i + 3ii'pi+r^pi = 0, (4.7a) 
P2 + SHp2 + rpQp2 = 0, (4.7b) 

PR + 4iJpR - r^pi - rpQP2 = 0, (4.7c) 

fiL + SHriL- SiTj^^Pi/mj^ = 0, (4.7d) 
+ SHriQ - {n'^'^f = 0. (4.7e) 

Here n^^ = C(3)T'^/7r^ is the equilibrium number density of the bosonic relativistic species; is a 
colUsion term for G production which, in the Umit of the massless gauginos, turns out to be [32, 33] 

CS- ^\ .E'^»?'°f-1 where Jfe) = PV^- 27, 72) 
° 16C(3)74fcr \9iJ = (1.634, 1.312, 1.271) 

and gi (with i = 1, 2, 3) are the gauge coupling constants of the MSSM calculated as functions of the 
temperature T. The latter quantity and the entropy density, s, can be obtained through the relations 

7^2 277-2 

PR = -9.T^ and s = —g,T\ (4.9) 
Also the Hubble expansion parameter, H, during this period is given by 

H = -j}—{mQnQ + pi+P2 + PRY^\ (4.10) 

Clearly, in the limit of massless MSSM gauginos, the computation is m^-independent. 

The numerical integration of Eqs. (4.7a)-(4.7e) is faciUtated by absorbing the dilution terms. To 
this end, we find it convenient to define [44] the following dimensionless variables 

/i = piR^ /2 = P2R^ /r = PrR\ h = riLR^ and /g = n^R^ (4.11) 

Converting the time derivatives to derivatives w.r.t. t = In with i?; being the value of the scale 

factor at the onset of the N oscillations - the precise value of Ri turns out to be numerically irrelevant 
- Eqs. (4.7a)-(4.7e) become 

Hf[ = -r^/i, (4.12a) 

Hf^ = -rpQ/2, (4.12b) 

Hfk = rj^hR + TpQf2R, (4.12c) 

Hfl = SlT^R^ (4.12d) 

HfL = C^irf'^fRK (4.12e) 

Also H and T can be expressed in terms of the variables in Eq. (4. 1 1) as 



JmQf^, + fi + f2 + fR/R / 30 



1/4 

H=- " , and r = I 7 ■'1^, ) • (4.13) 
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Figure 1: The evolution of the quantities log pi with i = 1 (dark gray line), i = 2 (gray line), i = K 
(light gray line), log p2i (dashed gray line), log Yl (black solid line) and log Yq (black dashed line) as 
functions of logT, for X = 0.0071, Aa = 10-^/a = lO^^ GeV, ii = 1 TeV, h^s = 10-^ kx = 1 
and On ^ 307, resulting to Fg = 2 • IQ-^^ andFz, = 2.5 • 10-^° formes = 0.05 eV and Jeff = 0.01. 



The system of Eqs. (4.12a)-(4.12e) can be solved, imposing the following initial conditions (the quan- 
tities below are considered functions of the independent variable t): 

PiiO) = pii, pr(0) = ng(0) = niiO) = and /92(Ti?PQ) = Pa, (4.14) 

where Xhpq is the value of t corresponding to the temperature Thpq which is defined as the solution of 
the equation H (Thpq) = mpq and can be found numerically. Needless to say that we set P2(t^) = 
for T < xhpq- 

In Fig. 1 , we illustrate the cosmological evolution of the quantities log pi with i = 1 (dark gray 
line), i = 2 (gray line), i = R (light gray line), log p2i (dashed gray line), log Yl (black soUd line) and 
log Yq (black dashed line) as functions of log T for values of the parameters which are allowed by all 
the restrictions described in Sec. 5. In particular we set A = 0.0071, Aa = 10~^, fa = 10^^ GeV, /x = 
1 TeV, hcs = 10"^ kx = I, cn ^ 307, = 0.05 cV and ^cff = 0.01. From Fig. l-(a), we 
observe that non-MiVl is followed by a matter dominated (MD) era, due to the oscillating and decaying 
inflaton system, which lasts until T c:iTf;j given by Eq. (4.2). The completion of the reheating process 
corresponds to the intersection of pi with pR. Afterwards the universe enters the conventional radiation 
dominated (RD) epoch of standard Big Bang cosmology. This becomes possible thanks to Eq. (4.6), 
since then the decay of the PQ system commences at T = Tj/pq and is completed at T ~ TpQ given 
by Eq. (4.5), before its domination over radiation. So a second episode of reheating does not occur. As 
we show below, this is a generic feature of our model. Due to this fact, the L and G yields, Yl = ul/s 
and = n^/s respectively, take their actual values (2.5 • lO"^'^ and 2 • 10~^^) for T ~ as shown 
in Fig. 2-(b). Both numerical values are compatible with observational data - see Sec. 5. 

4.3 Analytical Approach 

The numerical findings above can be understood by some simple analytic formulas. Most of them are 
widely employed in the literature - cf. Ref. [19, 22-25, 30]. In particular, the B yield can be computed 
as 



(a) Yb = -0.35Yl with (b) Yl = ul/s = clTj^ where cl = -5eL/4.m^. (4.15) 



5 Constraining the Model Parameters 
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The numerical factor in Eq. (4.15a) is due to sphaleron effects, whereas the unusual numerical factor 
in the definition of cl is due to the slightly different calculation [44] of - cf. Ref. [19]. On the other 
hand, the G yield at the onset of nucleosynthesis is estimated to be 

= n^/s ~ cgT^ with = 1.9 • IQ-'^^/GeV. (4.16) 

Both Eqs. (4.15) and (4.16) calculate the correct values of the B and G abundances provided that no 
entropy production occurs for T < T^. We show in the following that this is the case for our model. 

The evolution of the various energy densities involved in the post-inflationary dynamics can be 
well approximated - see e.g. Ref. [34,44] - by the expressions 

Pi = Piie-^^ for r>T^, (4.17a) 

P2 = Pae-^^^-^^P'^) for r>TpQ, (4.17b) 
PK = PK{T^){T/Tf^f for T<T^. (4.17c) 

Possible domination of p2 may occur for T < T^, since pR is steeper than p2. The equality between 
these two energy densities could be attained for T = Teq where 

PR(Teq) = P2(Teq) => T,^ ~ = (4.18) 

In deriving the equation above, we use the fact that pR(r^) = pi{T^) given by Eq. (4.17a) and that 
~ for Tcq <T< Tj^ due to the isentropic expansion. Also we assume that for T > Thpq we 
have a MD era driven by pi. This is a natural assumption since the condition of the out-of-equilibrium 
decay of N gives an upper bound on hes, which prevents the unlimited increase of T^. Indeed 



Pm^ _ 2\/2A7r 



" " CT V CT 

where Eqs. (4.2) and (4.1) are employed. The domination of p2 over the several energy densities - and 
therefore, a second reheating process - can be avoided if we impose the condition 

TpQ > Teq ^ A„ > \hlsM^c^/3V2Xl. (4.20) 



Combining the last relation with Eq. (3.10), we arrive at h^s < \/ 3\/2A^/MpqC7^. Taking e.g. 
\n ~ 0.01, cti ~ 10^ and MpQ ~ 10~^ the last relation implies heS ^ 100 which is much less 
restrictive than Eq. (4.19) which gives heS ^ 0.01 for A < 3.5 and ct ^ 0.3 - see Sec. 5.2. Therefore, 
the decay of the PQ system before its domination over radiation can be naturally accommodated within 
our set-up. Moreover, we remark that the results on T^, Yl and Yq are independent of phpq and p2i 
or yu, Aa and /„ when Eq. (4.20) holds. 

5 Constraining the Model Parameters 

We exhibit the constraints that we impose on our cosmological set-up in Sec. 5.1, and deUneate the 
allowed parameter space of our model in Sec. 5.2. 

5.1 Imposed Constraints 

Under the assumption that (i) the curvature perturbations generated by a are solely responsible for the 
observed curvature perturbations and (ii) the violation of Eq. (3.10) occurs after the violation of the 
slow -roll conditions in Eqs. (3.12a) and (3.12b), the parameters of our model can be restricted once we 
impose the following requirements: 
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5.1 .1 According to the inflationary paradigm, the horizon and flatness problems of the stan- 
dard Big Bag cosmology can be successfully resolved provided that A^* defined by Eq. (3.16) takes a 
certain value, which depends on the details of the cosmological scenario. Employing standard meth- 
ods [5,45], we can easily derive the required iV* for our model, consistently with the fact that the PQ 
oscillatory system remains subdominant during the post-inflationary era. Namely we obtain 

iV.^22.5 + 21n^^-^--ln^^ + -ln3-^ + -ln— (5.1) 



5.1 .2 The mflationary observables denved m Sec. 3.2 are to be consistent with the fittmg [27] 

1/2 

of the WMAP7, BAO and Hq data. As usually, we adopt the central value of , whereas we allow 
the remaining quantities to vary within the 95% confidence level (c.l.) ranges. Namely, 



(a) P^^ ~ 4.93- 10-^ (b) ris = 0.968 ±0.024, (c) -0.062 < < 0.018 and (d) r < 0.24. (5.2) 

5.1 .3 For the realization of non-MTVI , we assume that c-jz takes relatively large values - see 
e.g. Eq. (3.5). This assumption may [4, 46] jeopardize the validity of the classical approximation, 
on which the analysis of the inflationary behavior is based. To avoid this inconsistency - which is 
rather questionable [8, 46] though - we have to check the hierarchy between the ultraviolet cut-off, 
A = mp/cTz, of the effective theory and the inflationary scale, which is represented by V^j(c*)i/^ or, 
less restrictively, by the corresponding Hubble parameter, if* = Vj;f^{a^,)^^'^ /VSmp. In particular, the 
validity of the effective theory implies [46] 

(a) %j(cj*)i/^ < A or (b) < A for (c) c-jz > 1- (5.3) 

5.1.4 In agreement with our assumption about hierarchical light neutrino masses and the 
results of neutrino oscillation experiments [47], m^s - involved in the definition of in Eq. (4.4) - 
can be related to the squared mass difference measured in atmospheric neutrino oscillations, Am^. 
Taking the central value of the latter quantity, we set 



m^3 ~ y^Am| = (2.43 • lO"''^) ^''^ cV ~ 0.05 eV. (5.4) 

This value is low enough to ensure that the lepton asymmetry is not erased by lepton number violating 
2^2 scatterings [48] at all temperatures between and 100 GeV. 

5.1 .5 The interpretation of BAU through non-thermal leptogenesis dictates [27] at 95% c.l. 

Yb = (8.74 ± 0.42) • 10-11 => 8.32 < Yb/IO-h < 9.16. (5.5) 

Given our ignorance about ^eflf in Eq. (4.4), we impose only the lower bound of the inequality above 
as an absolute constraint. 

5.1.6 In order to avoid spoiling the success of the SBB nucleosynthesis, an upper bound 
on Yq is to be imposed depending on the G mass, mg, and the dominant G decay mode. For the 
conservative case where G decays with a tiny hadronic branching ratio, we have [33] 



'lo- 


-15 


10- 
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-14 


10- 


-13 


10- 

V 


-12 



for rriQ 2± 



0.45 TcV 
0.69 TeV 

(5.6) 

10.6 TeV 
13.5 TeV. 



The bound above can be somehow relaxed in the case of a stable G. However, it is achievable in our 
model, as we see below. 
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(a) (b) 

Figure 2: The allowed by Eqs. (5.1), (5.2a), (5.3b) and (5.3c) values of A (solid line) and T^h - 
given by Eq. (4.2) - (dashed line) [ai (solid line) and cr* (dashed line)] versus cn (a) [(b)] for Aa = 
10~^, fcx = 1 andhes = 10~^. Also,a2c (dot-dashed line) given by Eq. (3.10d) as function of ctz is 
shown. The light gray and gray segments denote values of the various quantities satisfying Eq. (5.3a) 
too, whereas along the light gray segments we obtain cr* > mp. 



5.2 Results 

As can be easily seen from the relevant expressions above, our cosmological set-up depends on the 
following parameters: 

A, Aa, fa, A^, kx, c-jz, hes and ^efr- 

Our results are independent of Aa and kx, provided that Eqs. (3.10) and (4.20) are satisfied. With 
these conditions, the contribution of T4c to V^j remains subdominant and the PQ system decays before 
radiation domination. We therefore set A^ = 10~^ and kx = 1 throughout our calculation. The chosen 
Aa is close to its largest value allowed by Eq. (3.10b), whereas kx is fixed to a natural value. Also ^cff 
affects exclusively the Yl caclulation through Eqs. (4.4) and (4.15). We take S^s = 1, which allows 
us to obtain via Eq. (4.4) the maximal [49] possible lepton asymmetry. This choice in conjuction with 
the imposition of the lower bound on Yg in Eq. (5.5) provides the most conservative restriction on our 
parameters. Also, we set A^ = 0.01 [A^ = 1] so as to obtain /i ~ 1 TeV with = 10^^ GeV [/„ = 
10^^ GeV] - evidently, the generation of the fi term of the MSSM through the PQ symmetry breaking 
does not favor lower values for fa. As we show below, the selected values for the above quantities give 
us a wide and natural allowed region for the remaining fundamental parameters (A, c-ji and hes) of 
our model. In our numerical code, we use as input parameters cr*, hcs, fa and ct^. For every chosen 
ct^ > 1 and h^s, we restrict A and a* so as the conditions Eq. (5.1) - with T^h evaluated consistently 
using Eq. (4.2) - and (5.2a) are satisfied. Let us also remark that in our numerical calculations, we use 

1/2 

the complete formulas for the slow-roll parameters and P-^' in Eqs. (3.12a), (3.12b) and (3.18) and 
not the approximate relations, which are listed in Sec. 3.2 for the sake of presentation. 

Our results are presented in Fig. 2, where we draw the allowed values of A (solid line) and Tj-h 
(dashed line) [cxf (sohd line) and a* (dashed line)] versus c-ji (a) [(b)] for heS = 10~^ and fa = 
10^2 GeV. In Fig. 2-(b) we also draw a2c - derived from Eq. (3.10b) - as a function of ctz. The upper 
[lower] bound on c-ji comes from the saturation of the inequality in Eq. (5.3b) [Eq. (5.3c)]. On the 
other hand, Eq. (5.3a) is vahd along the gray and hght gray segments of the various curves. Along the 
light gray segments, though, we obtain cj* > mp. The latter regions of parameter space, although can 
be considered as less favored, are not necessarily excluded [26], since the energy density of the inflaton 
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remains sub-Planckian and so, corrections from quantum gravity can be assumed to be small. In all, 
we obtain 

l<cn<l.56-10^ and 2.6 ■ IQ-^ < A < 3.5 for 52.5 < iV, < 54.6. (5.7) 

From Fig. 2-(a), we observe that A depends on c-n almost linearly. This can be understood by combining 
Eq. (3.18) and Eq. (5.2a). The resulting relations reveal that A is to be proportional to c-jz, so as 
Eq. (5.2a) is satisfied with almost constant A^*. Indeed we find 

A = 3 • 10-\cTz/N^ =^ cn = 41850A for iV* ~ 55. (5.8) 

On the other hand, the variation of Cf and a^: as a function of ctz - drawn in Fig. 2-(b) - is consistent 
with Eqs. (3.17) and (3.14). If ctz varies within its allowed region as given in Eq. (5.7), we obtain 

0.963 <n,< 0.965, -6.8 < ^ < -6.1 and 4.4 > ^ > 3.4. (5.9) 

Clearly, the predicted ng and r lie within the allowed ranges given in Eq. (5.2b) and Eq. (5.2c) respec- 
tively, whereas ccg remains quite small. These findings depend very weakly on Tj-h - and therefore on 
hes, since this controls the value of Tj-h via Eq. (4.2) - because Tj-h appears in Eq. (5.1) through the 
one third of its logarithm, and consequently its variation upon some orders of magnitude has a minor 
impact on the required value of A^* . 

On the contrary, plays a key-role in simultaneously satisfying Eqs. (4.19), (5.5) and (5.6) - see 
Eqs. (4.15) and (4.16). For this reason we display in Fig. 3-(a) [Fig. 3-(b)] the allowed regions by all 
imposed constraints in the A - h^s plane for fa = 10^^ GeV [fa = 10^^ GeV] - cf. Ref. [23]. The 
restrictions on the parameters arising from the inflationary epoch are denoted by dotted and double- 
dotted dashed lines, whereas the ones originating from the post-inflationary era are depicted by solid, 
dashed and dot-dashed lines . In particular, the double-dotted dashed [dotted] lines come from the 
bounds of Eq. (5.3b) [Eq. (5.3c)]. In the horizontally fined regions Eq. (5.3a) holds, whereas in the 
vertically hatched region we get cr* > mp. On the other hand, the solid [dashed] lines correspond to the 
lower [most conservative upper] bound on Yb [Ig] in Eq. (5.5) [Eq. (5.6)]. Since we use \Ses\ = 1, it is 
clear from Eqs. (4.4) and (4.15) that values of heS above the solid fine are compatible with the current 
data in Eq. (5.5) for conveniently adjusing \Ses\ < 1- However the strength of hcs can be restricted 
by the bounds of Eq. (5.6), which can be translated into bounds on via Eq. (4.16). Specifically we 
obtain ~ (0.1 - 1) • lO'^^ or T^h ^ (0.53 - 5.3) • 10^ GeV (gray area), Y^ ~ (0.1 - 1) • lO'^^ 
or Trh (0.53 - 5.3) • 10^ GeV (yellow area), Y^ ~ (0.1 - 1) • lO-^'' or T^h (0.53 - 5.3) ■ 10^ GeV 
(dark gray area) and Yq ~ (0.1 - 1) • 10"^^ or T^h c± (0.53 - 5.3) • 10^ GeV (fight gray area). The 
competition of the two restrictions above can be presented also analytically. Indeed, plugging Trh via 
Eq. (4.2) into Eqs. (4.15) and (4.16), seting Yl = y/""™ and = yif^^ - the exact numerical values 
can be extraced by Eqs. (5.5) and (5.6) respectively - and solving the resulting equations w.r.t. h^s, we 
obtain the following inequifities: 



max 



2 /tt ymin 2 Y~ 

(a) > , ^ ^V- and (b) < , ^ (5.10) 

These reproduce quit accuratelly the behavior seen in Fig. 3. On the other hand, the out-of-equifibrum 
condition depicted by a dashed fine - see Eq. (4.19) - puts the upper bound on /leff in a minor portion 

of the parameter space. 

Comparing Fig. 3-(a) and Fig. 3-(b), we conclude that given by Eq. (4.1a) is reduced for 
fa = 10^^ GeV w.r.t. its value for /„ = 10^^ GeV, and so the condition of Eq. (4.19) cuts a larger slide 
of the available parameter space. Letting A vary within its allowed region in Eq. (5.7), we obtain 

4.3 • 10"^ [1.3 • 10-^] < /ieff/10"^ < 0.58 for /„ = 10^^ GeV [/« = 10^^ GeV] (5.11) 
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N a 

Figure 3: Allowed (shaded) regions as determined by Eqs. (4.19), (5.3b), (5.3c), (5.5) and (5.6) in 
the A - /loff plane, for Aa = 10~^ A;^ = 1, M = 1 TeV and fa = lO^^ GeV (a) or fa = 10" GeV 
(b). The conventions adopted for the various lines and shaded or hatched regions are also shown. 



where the overall minimal [maximal] /icfr can be found in the upper, almost central [lower right] corner 
of the allowed region. As we see above and can be induced by Eqs. (4.19) and (5.10b), the maximum 
allowed h^s is fa independent and it is obtained for A ~ 0.005 [A ~ 0.05] and fa = 10^^ GeV [fa = 
10^^ GeV]. This point gives also a lower bound on \6cs\, |5eff| ^ 2 • 10^^ - which is obviously also fa 
independent. Note finally that in both cases the resulting m^'s can be much lower than those obtained 
within the simplest model of sneutrino inflation with a quadratic potential [21,22]. On the other hand, 
/leff turns out to be comparable with the one obtained in those models. 

6 Conclusions 

In this paper we attempted to embed within a realistic cosmological setting one of the recently formu- 
lated [8] SUSY models of chaotic inflation with non-minimal coupling to gravity. We concentrated 
on a moderate extension of the MSSM augmented by three RH neutrino superfields and three other 
singlet superfields, which lead to a PQPT tied to renormalizable superpotential terms. The coupling 
between the RH neutrinos and one of the fields associated with the PQPT plays a crucial role for the 
implementation of our scenario. We showed that the model non only supports non-MI driven by the 
lightest RH sneutrino, but it also resolves the strong CP and the /i problems of the MSSM and, even 
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more, it leads to the production of the required by the observations BAU via non-thermal leptogenesis, 
which accompanies the inflaton's decay. Moreover the G abundance becomes observationally safe for 
G masses even lower than 10 TeV. An important prerequisite for all these is that the parameter of 
the superpotential related to the PQPT, Aa, is adequately small. Imposing a number of observational 
constraints arising from the data on the inflationary observables, the BAU, the concentration of the 
unstable G at the onset of nucleosynthesis and the mass of the heaviest light neutrino, we restrict the 
effective Yukawa coupling, involved in the decay of the inflaton, to relatively small values, and the 
inflaton mass to values lower than 10^^ GeV. 
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Appendix: Non-Minimally Curvature-Coupled Scalars in SUGRA 

Non-MI can be realized by a scalar field with a non-minimal coupling to the Ricci scalar curvature. The 
formulation of a such theory within SUGRA is described below. Recall that we follow the conventions 
of Ref. [34] for the quantities related to the gravitational sector of our set-up. 

In contrast to the non-SUSY case - see e.g. Ref. [1,2,5,9] - we find it convenient to start our analy- 
sis with the general Einstein-frame (EF) action for the scalar fields plus gravity in four dimensional, 
A/" = 1 SUGRA [7,8]: 



/ 



d^x^ {-\mln + K^pgi'^D^<i)'^D,<t>*~^ - , (A.l) 



where hat is used to denote quantities defined in the EF; g is the determinant of the Friedmann- 
Robertson- Walker background metric [34]; 

""-'^ = > ° = ^'^^^ " ^^^^ ^^-^^ 

are the covariant derivatives for scalar fields Here stand for the vector gauge fields and fe^ is 
the Killing vector, defining the gauge transformations of the scalars [8]. Assuming that the D-terms 
of (/)" vanish - as for the singlet scalars N,P,X and X in our model - the EF scalar potential, V, is 
given in terms of the Kahler potential, K, and the superpotential, W, by 

V = e^/^^p (^K^'^FaF^ - ' ^^-^^ 

with K'^^K^/s = 5^ and = W^^^ + K^^cW/m^,. 

The action in Eq. (A.l) can be brought the Jordan frame (IF) by performing a conformal transfor- 
mation [38]. Indeed, if we define the JF metric, gf^l,, through the relation 

- ^ ('v^=J^2V=5/9 and ?''- = -35''7J2, 

1 /2 

- where □ = {-g)~ ' iV^d'') 

- we obtain the action in the JF as follows 
S = j d^x^ (^^^ + it^M^^"^ - hlK^-^D^^^D^'<t>*~^ - with V = ^V. (A.5) 
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Taking into account that d^Q = D^^Q - since Q is gauge invariant - and that the purely bosonic part, 
Afj,, of the on-shell value of the auxiliary field is given by 

= -^rnl {D^cl^'^na - D^r^na) (A.6) 
- with JIq, = ii^^a and il^ = ^,</>*« - and for the choice 

n = -3e-^/^"'p ^ K = -3ml In (-^^/3) , (A.7) 
we arrive at the following action 

S = j (^^^ + 'rnl^^-pDf,(j)'^Di'^*~^ - ^A^.Ai'/ml - . (A.8) 

It is clear from the first term of the RHS of this expression that the resulting S exhibits non-minimal 
couplings of the (/)"'s to IZ. However, $7 enters in the kinetic terms of the (/)"'s too. In order to get 
canonical kinetic terms, we need Q.^^ = 6^^ and A^ = 0. The first condition is satisfied [8] by the 
choice 

n = -3 + 5^p^^ - 3{¥{r) + F*(0*^)), (A.9) 

where F is a dimensionless, holomorphic function, which expresses the non-minimal coupling to grav- 
ity. Note that even when F(^") = for some a, the 0"'s are conformally coupled to gravity due to the 
second term of the RHS of the expression above. This choice for the frame function leads via Eq. (A.7) 
to the following Kahler potential 

K = -3ml In (l - ^'^a^'/'^*'^ + + F* (</>*"))) • (A.IO) 

On the other hand, A^ = when the dynamics of the ^"'s is dominated only by the real moduli 
|<^"|. Therefore, it is possible to get a SUGRA realization of the inflationary models with non-minimal 
coupling. 
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